Mon. Not. R. Astron. Soc. 000,[T]{T5](2010) Printed 22 January 2010 (MN MTeX style file v2.2) 



Keck telescope constraint on cosmological variation of the 
proton-to-electron mass ratio 



o 



A. L. Malec/* R. Buning,^ M. T. Murphy/ N. Milutinovic,^ S. L. Ellison,^ 
J. X. Prochaska,"^ L. Kaper,^'^ J. Tumlinson,^ R. F. Carswell/ W. Ubachs^ 

^ Centre for Astrophysics and Supercomputing, Swinburne University of Technology. Melbourne, Victoria 3122, Australia 

^Laser Centre, VU University, De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 

^ Department of Physics and Astronomy, University of Victoria, Victoria, BC, V8P lAl, Canada 

* University of California Observatories - Lick Observatory, University of California, Santa Cruz, CA 95064 

^Astronomical Institute Anton Pannekoek, Universiteit van Amsterdam, 1098 SJ Amsterdam, The Netherlands 

^Yale Center for Astronomy and Astrophysics, Department of Physics, New Haven, CT 06520, USA 

''institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 OHA, UK 



o 
u 

6 



> 

oo 

o 

O 
O 



Accepted 2009 December 14. Received 2009 November 18; in original form 2009 June 9 



ABSTRACT 

Molecular transitions recently discovered at redshift Zabs - 2.059 toward the bright back- 
ground quasar J2123-0050 are analysed to limit cosmological variation in the proton-to- 
electron mass ratio, jj. = nip/me. Observed with the Keck telescope, the optical echelle spec- 
trum has the highest resolving power and largest number (86) of H2 transitions in such anal- 
yses so far Also, (seven) HD transitions are used for the first time to constrain /^-variation. 
These factors, and an analysis employing the fewest possible free parameters, strongly con- 
strain yu's relative deviation from the current laboratory value: Ayu/ju - (+5.6±5.5stat±2.9sys)x 
10"^, indicating an insignificantly larger ju in the absorber This is the first Keck result to com- 
plement recent null constraints from three systems at Zabs > 2.5 observed with the Very Large 
Telescope. The main possible systematic errors stem from wavelength calibration uncertain- 
ties. In particular, distortions in the wavelength solution on echelle order scales are estimated 
to contribute approximately half the total systematic error component, but our estimate is 
model dependent and may therefore under or overestimate the real effect, if present. 

To assist future //-variation analyses of this kind, and other astrophysical studies of H2 in 
general, we provide a compilation of the most precise laboratory wavelengths and calculated 
parameters important for absorption-line work with H2 transitions redwards of the hydrogen 
Lyman limit. 
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1 INTRODUCTION 

The Standard Model of particle physics is parametrized by several 
dimensionless 'fundamental constants', such as coupling constants 
and mass ratios, whose values are not predicted by the Model it- 
self. Instead their values and, indeed, their constancy must be estab- 
lished experimentally. If found to vary in time or space, understand- 
ing their dynamics may require a more fundamental theory, perhaps 
one unifying the four known physical interactions. Two such pa- 
rameters whose constancy can be tested to high precision are the 
fine-structure constant, a = e^/Aneohc, characterising electromag- 
netism's strength, and the proton-to-electron mass ratio, yj = m^/m^ 
- effectively the ratio of the strong and electro- weak scales. Grand 
Unified Theories can predict relationships between variations in /j 
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and a (e.g. ICalmet & Fritzscl]|2002l) . but which parameter varies 
the most and whe ther in the same o r opposite sense to the other is 
model-dependent l lDent et alJlOOj) . 

Earth-bound laboratory experiments, conducted over several- 
year time-scales, which use ultra-stable lasers to compare dif- 
ferent atomic clocks based on different a toms/ions (e.g. Cs, 



Hg+ A1+, Yb+ Sr, Dy; e.g. |P restage et alj 1 19951 ; iMarion et all 
I2OO3I ; IPeik etai]|2004l ; Icingbz et al. 2007), have limited a's time- 
deriva tive to a/a = (-1.6 ± 2.3) X 10"'^ yr"' l lRosenband et al.l 
I2OO8I) . Combining results from similar experi ments limits fi's ti me- 
derivative to ft/fi = (-1.6 ± 1.7) x 10 '^^ yr"' jBlatt et alj|2008l) . A 
direct constraint was also recently obtained from a ro-vibrational 
transition in the SFe mol ecule, fi/fi = (-3.8 ± 5.6) x 10"''' yr"' 
jShelkovnikov et alj|2008h . 

Important probes of variations over much larger space- and 
time-scales - up to ~90% of the age of the Universe - are nar- 
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row absorption lines imprinted on the spectra of distant, back- 
grou nd quasars by gas c louds associated with foreground galax- 
ies teahcall et alj Il967h . In particular, the Lyman and Werner 
transitions of molecular hydrogen (laboratory wavelengths /ii^i, < 
1150 A), by far the most abundant molecule in the Universe, are 
useful ^-variation indicators at cosmological red shifts z> 2 where 
they are detectable with ground-b ased telescopes jThompsonl 1 9751 : 
IVarshalovich & Levshakovlll993h . 

Variations in ji should shift the ro-vibronic transition frequen- 
cies in molecular spectra. This mass-dependent shift is quantified 
by a sensitivity coefficient. A",, for each transition /. Consider a 
single transition arising in an absorption cloud whose redshift is 
established to be z^i,, from other transitions which are insensitive 
to variations in ji. If /j was the same in the absorption cloud as 
in the laboratory, we would expect to find transition ; at wave- 
length X^^(l + Zabs)- If instead we measure it to be at wavelength 
Ai = + Zi) (i.e. at redshift z, Zabs). then the shift in redshift, 

Az, = Zi - Zab.s, or velocity, Av,, can be ascribed to a variation in yu, 

^.-^ = K,^, (1) 

c 1 + Zabs y" 

where Ayu//j = (jj.- /Jiab)/y"iab for y"iab and yU- the current laboratory 
value of fi and its value in the absorption cloud at redshift z, re- 
spectively. That different transitions have different K values enables 
a differential measurement of Aji/fi from two or more transitions. 
That is, it allows the redshift Zabs to be determined simultaneously 
with Ayu//j. 

Indications for a significantly positive Ayu/yU have been derived 
from two H2-bearing quasar absorbers observed with the Ultravio- 
let and Visual Echelle Spectrograph ( UVES) on the ESQ K ueyen 
Very Large Telescope (VLT) in Chile. Ilvanchik et al] ( |2005|) stud- 
ied high resolving power (R ^ 53000), high signal-to-noise ratio 
(SNR « 30-70) UVES spectra of the Zab.s = 2.595 and 3.025 
absorbers towards the quasars Q 0405-443 and Q 0347-383, re- 
spectively. They used a total of 76 Ht lines in the two spec- 
tra, fitting them 'line-by-line' - i.e. independently of each other 
- to derive values for /I,. In this analysis the H2 absorption pro- 
files in Q 0347-383 were treated as comprising a single cloud. Of 
the two resolved H2 features observed in Q 0405-443, only the 
stronges t was fitted. Using tw o di fferent sets o f Hi laborat ory wave- 
lengths jAbgrall etal]|l993l and IPhilip et~al] 12004 ) Ivanchik et al.l 
derived two sets of z, values which yielded two Aji/fi values, 
{+3 5 ± 7.5) X 10'^ and ( +16.5 ± 7.4) x 10^*^ respectively. 

iReinhold et alj ^2W^l) subsequently improved the laboratory 
wavelengths and the calculation of the K sensitivity coefficients 
(see Sections 12.21 and 123). Using the same values of /I, and t heir 
uncertainties derived by Ilvanchik et alJ ( l2005h . lReinhoId et al.| per- 
formed a 'line-by-line' analysis to refine the values of Ayu/yU for 
each absorber: (+20.6 ± 7.9) x 10^* for Q 0347-383 and (+27.8 ± 
8.8) X lO"*- for Q 0405-443. The combined value of (+24.5 ± 5.9) x 
10"* was presented as an indication for cosmological variation in 

iKing et al.l ( I2OO8I) reanalysed the same raw UVES quasar 
spectra with improved flux extraction and, more importantly, us- 
ing the improved w a velength cal ibration procedures detailed in 
iMurphv et all ilOOH) . iKing et al.l used slightly more H2 transi- 
tions and a 'simultaneous fitting' technique with explicit treatment 
of Lyman-Q' forest lines (see Section 13. Il l to constrain the val- 
ues of A/i/yU in the newly reduced and calibrated spectra. They 
found decreased values in both absorbers compared to previous 
works: Ayu/// = (+8.2 ± 7.4) x 10"' and (+ 10.1 ± 6.2) x 10"'^ for 
Q 0347-383 and Q 0405-443 respectively. iKing et alj also anal- 



ysed 64 H2 lines of a third absorption system, that at Zabs = 2.811 
towards Q 0528-250, using UVES spectra with R « 45000 and 
SNR a; 25^5. This provided the tightest constraint of all three 
absorbers: Ayu/^ = (-1.4 ± 3.9) x 10"'. Thus, the combined result, 
where the slightly positive values for Q 0347-383 and Q 0405-443 
are somewhat cancelled by the slightly negative but more precise 
value for Q 0528-250, was a null constraint of A/j/yu = (+2.6 ± 
3. 0) X 10"'. The same UV ES s pectra were also recen tly studied 
bv lWendt & ReimersI |2008') and lThompson et alj ( [20091) using dif- 
ferent data reduction and analysis techniques, generally aimed at 
avoiding and/or understanding potential systematic errors and bi- 
ases in Ayu/yu. They also fin d null constr aints, albeit with somewhat 
larger statistical errors than lKing et al.l due to their more conserva- 
tive approaches. 

So, whil e som e analyses indicate a varying yu (e.g. 
IReinhold et all l2006h , recent reanalyses of the three well- 
documented and high-quality H2 absorption spectra at z > 2 do 
not provide evidence for cosmological variation in yu. A much 
larger statistical sample is desirable, deriving from several tele- 
scopes and spectrographs, for definitive results and to provide mea- 
surements over a larger redshift range. However, the scarcity of 
known Ht absorbers hampers such progress: >1000 absorption sys- 
tems rich in neutral hydrogen - i.e. damped Lyman-O' systems, 
with Hi column densities A^( Hi) > 2 x 10^" cm"^ - are known 
(e.g. lProchaska & Wolfel2o6i) but syste matic searches for H2 have 
been conducted in <I00 systems (e.g. iLedoux et al] l2003h and 
only ~15 are known to h arbour detectable column densities of Ht 
jNoterdaeme et al ]|2008l) . Furthermore, for varying-yu analyses, Ht 
absorbers must be (i) at z > 2 to shift enough Lyman and Werner 
transitions above the atmospheric cutoff (~3000 A), (ii) have bright 
background quasars to enable high-SNR, high resolution spec- 
troscopy and, (iii) have high enough H2 column densities so that 
individual transitions absorb significant fractions of the quasar con- 
tinuum. 

In these respects, and others, the newly discovered absorber 
studied here - at Zabs = 2.059 towards the Zcm = 2.261 quasar 
SDSS J2 12329.46-005052.9 (hereafter J2 123-0050; Milutinovic 
et al., in preparation) - is exceptional. It is the firs t studied with the 
High Resolution Echelle Spectrometer (HIRES: pN^gt et alj[T994l) 
on the Keck telescope in Hawaii to provide constraints of similar 
precision to the three VLT onefl J2123-0050 is unusually bright 
(r-band magnitude a 16.5 mag), providing a high SNR spectrum 
in just 6 hours observation with, importantly, the highest spectral 
resolution for an H2 absorber to date: R 1 10000 or full-width- 
at-half-maximum FWHM a 2.7 km s"'. HIRES's high ultravio- 
let (UV) throughput provides 86 H2 transitions for constraining 
yu-variation, the largest number in an individual absorber so far. 
Also, for the first time, (7) HD transitions are used to constrain 
^-variatio n; HP has been observed in just two other high-re dshift 
absorbers jVarshalovich et al.ll200ll : lNoterdaeme et al.F 2008Y 

This paper is organised as follows. The data required for our 
analysis - the Keck spectrum of J2123-0050, laboratory H2/HD 
transition wavelengths and the sensitivity coefficients - are de- 
scribed in Section|2] Section[3]presents our analysis technique and 
main result while in Section|4]we conduct a variety of internal con- 
sistency checks and consider the most important systematic errors. 



' ICowie & Songailal i 19951) previously analysed Keck/HIRES spectra of 
Q0528-250 to derive a comparatively weak constraint of Afi/jj = (-80 ± 
313) x 10"'. 
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We conclude with a comparison of our results with others in the 
literature in Section[5] 



2 DATA 

2.1 Keck spectrum of J2123-0050 

Full details of the Keck/HIRES observations of J2 123-0050 are 
provided in Milutinovic et al. (in preparation) who study the phys- 
ical conditions and elemental abundances in the z-jos = 2.059 ab- 
sorber. Here we present only the features of the spectral data reduc- 
tion and the final spectrum which are important for the varying-// 
analysis. 

The HIRES observations comprised 2 sets of 3x1 -hr expo- 
sures of J2123-0050 taken on consecutive nights (2007 August 19 
& 20). The seeing was very good, just 0'.'3-0'.'5 at ~ 6000 A (where 
the guide camera is sensitive), and so a O'.'4-wide slit was used. A 
single thorium-argon (ThAr) lamp exposure calibrated each set of 
three quasar exposures; the ThAr exposure was taken immediately 
before or after each quasar set. Spectrograph temperature and at- 
mospheric pressure shifts between the quasar and ThAr exposures 
were < 1 K and < 1 mbar respectively. 

The raw spectra were reduced using the hiredux softwar^l 
written and maintained by one of us (JXP). The flux extraction pro- 
cedure also calculates the formal statistical flux error spectrum for 
each echelle order. Particular attention was paid to accurate wave- 
length calibration. ThAr lines w e re pre -selected with the proce- 
dures described in iMurphv et aO ( l2007h and improvements were 
made to the ThAr line fitting in hiredux to ensure that reliable cen- 
troids for the ThAr lines were determined. The ThAr flux was ex- 
tracted using the same spatial profile weights as the correspond- 
ing quasar extraction. The root-mean-square (RMS) of the wave- 
length calibration residuals was ~80ms"'; that is, the calibration 
at any given wavelength has RMS error ~80ms"'. The resolv- 
ing power was measured from the extracted ThAr spectra to be 
R s! 110000. During extraction, the spectra were re-binned to a 
common vacuum-heliocentric wavelength scale with a dispersion 
of 1.3 km s"' pixer'. They were combined to form a final spec- 
trum for subsequent analysis using the uves_popler softwar^l writ- 
ten and maintained by one of us (MTM). 

The final spectrum of J2 123 -0050 covers (vacuum- 
heliocentric) wavelengths 3071-5896 A. An initial global quasar 
continuum was constructed using simple polynomial fits. The de- 
tails of this process are unimportant because the local continuum 
was modelled and fitted simultaneously with most H2/HD and 
Lyman-a forest absorption lines. That is, the initial global contin- 
uum is just a nominal, fixed starting guess for most molecular tran- 
sitions. All the H2/HD transitions fall bluewards of 3421 A where 
the SNR in the nominal continuum ranged from 7 per 1.3 km s"' 
pixel at ~3075 A to 25 at ~3420 A. We made a final, manual check 
on the data quality in the spectral region around each H2/HD tran- 
sition. In particular, we ensured that the data contributed by the dif- 
ferent exposures agreed to within the uncertainty expected from the 
individual flux error arrays derived during the extraction procedure. 
This is also a check that the statistical uncertainties fairly estimate 
the RMS flux variations in the final spectrum. Also, no significant 
wavelength or velocity shifts between the individual spectra were 
apparent. 

~ http://www.ucolick.org/~xavier/HIRedux 

' http://astronomy.swin.edu.au/~mmurphy/UVES_popler.html 
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Velocity around z = 2.0594 [km s 

Figure 2. Some of the 86 H2 and 7 HD lines from the J2123-0050 Keck 
spectrum on a velocity scale centred at z^bs = 2.0594. The spectrum (black 
histogram) is normalized by a nominal continuum (upper dotted line). Lo- 
cal linear continua (upper dashed lines) and zero levels (lower dashed lines) 
are fitted simultaneously with the molecular and broader Lyman-a lines, 
the positions of which are indicated with, respectively, lighter and darker 
tick marks (offset vertically from each other) above the spectrum. The 4- 
component fiducial fit is the solid curve. Note that the two left-most com- 
ponents are nearly coincident in velocity. 



Figure [T] shows the region of the final spectrum containing 
H2/HD lines (3071-3421 A) and a region redwards of the Lyman- 
a emission line (4905^926 A) containing metal lines which are 
required in our fit described in Section [3T| 

Figure|2]highlights some of the H2 and HD transitions cover- 
ing the observed range of overall line-strengths, SNRs and ground- 
state rotational levels, characterised by the quantum number J (V- 
levels'). H2 lines are observed in the B'E^-X'E;!^ Lyman and C'llu- 
X' Werner bands, for J e [0, 5] . HD is observed in six RO Lyman 
lines and one RO Werner line. The molecular absorption shows two 
distinct spectral features (SFs), separated by a 20kms"'. Figure[T] 
shows that the strong left-hand SF appears saturated for most low- 
J transitions, while the weaker right-hand SF appears unsaturated 
in almost all transitions. Only the left-hand SF is detected in the 
HD transitions. Clearly, accurately recovering A/j/yU with the maxi- 
mum precision available from the spectrum must take into account 
this non-trivial absorption profile, the strength of which varies from 
transition to transition, while also allowing the varying SNRs, de- 
grees of Lyman-Q- blending, and possible continuum-placement er- 
rors to contribute to the final uncertainty on Afi/ju. We achieve this 
with a simple x' fitting technique combined with a detailed model 
of the absorption profile in Section [3T| 

2.2 Laboratory wavelengths for H2 and HD transitions 

Equation (TJ states that velocity shifts between different H2/HD 
transitions can be related to a varying yu. But it assumes that the 
current laboratory wavelengths of those transitions are known to 
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Wavelength [A] 

Figure 1. All regions of the J2123-0050 Keck spectrum fitted simultaneously in our analysis. Here we show only 5 of the 37 panels comprising the full 
figure; the other panels are available in the electronic version of this paper. The spectrum (black histogram) is normalized by a nominal continuum (upper 
dotted fine) fitted over large spectral scales. Local linear continua (upper dashed lines) and zero levels (lower dashed lines) are fitted simultaneously with the 
H2/HD and broader Lyman-n lines. The fits are shown with solid grey/green lines. H2/HD transitions are labelled and their constituent velocity components 
are indicated by grey/green tick-marks immediately above the spectrum. Higher above the spectrum are tick-marks indicating the positions of Lyman-o- lines 
(dark-grey/blue) and Fen lines (lighter grey/red). Note that the metal-line velocity structure is constrained with the Fe nil 608 A transition shown in the final 
panel of the figure (i.e. in the electronic version). The residual spectrum (i.e. [data] - [fit]), normalized to the l-o" errors (faint, horizontal solid lines), is shown 
above the tick-marks. 
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high enough precision to be considered exact compared to the pre- 
cision available from the astronomical spectra. And with wave- 
lengths /liab £ 1150 A, laboratory measurements of these transi- 
tions are challenging. Nevertheless, the laboratory measurements 
have improved over the years and now provide accurate and pre- 
cise laboratory wavelengths for all the H2/HD transitions detected 
in our spectrum of J2 123-0050 and, more generally, for those rel- 
evant in /j-variation studies. 

While the classical data of the Meudon group reached frac - 
tional wavelength accuracies ofSAIA ~W-^ l lAbgrall etal]|l993h . 
the first laser calib ration study using a pulsed dye laser system 
( lHinneneta"i]|l994) reached similar or somewhat improved accu- 
racies of 6x 10"'. The implementation of Fourier tran sform-limited 
laser pulses in the harmonic conversion processes jUbachs et al.l 
Il997h allowed an order of magnitude improvement. Using this 
technique, several subsequent studies provided a data set of Ly- 
man and Werner transition wavelengths with fractional accuracies 
of 5x lO""* (Philip et al. 2004; Ubachs & Reinhold 200j). Recently, 
even more accurate values - fractional accuracies of 5 x 10"' for 
most Lyman transitions and 1-2 x 10"** for Werner transitions - 
were obtained from combining two-photon excitation results with 
Fourier-transform studies of emission lines between H2 excited 
states ISalumbides et al. 20oi). 

In our analysis we use the most accurate wavelength from 
these studies for each H2 transition detected in our spectrum of 
J2123-0050. Table [T] compiles the most precise laboratory data 
into a complete, up-to-date catalogue for all allowed Lyman and 
Werner H2 transitions between the lowest 8 rotational levels in 
the ground and excited states with wavelengths generally above 
the hydrogen Lyman limit (more specifically, to the first 20 and 
6 excited vibrational levels for Lyman and Werner transitions, re- 
spectively). Most important to /j- variation studies like ours are the 
wavelengths and sensitivity coefficients (see below), but we have 
also included the other laboratory data required for fitting the H2 
absorption lines seen in astronomical spectra, such as oscillator 
strengths. And while we detect only a small subset of these H2 tran- 
sitions in our spectrum of J2 123-0050, we include a complete list 
in Table [T] as a useful reference catalogue for other astronomical 
studies. 

The HD transitions were recently measured via direct ex- 
treme UV laser excitati on to a relative accuracy of a 5 x 10"** by 
iHollenstein et al. Illooi) and llvanov et al. and we use those 

values here. Table |2] provides the laboratory data for the 10 HD 
transitions falling in our spectrum, 7 of which we use in our analy- 
sis. 

2.3 Sensitivity coefficients, K 

The K co efficients used here were calculated via a semi-empirical 
analysis dUbachs et ai]|2007h . In this approach a Dunham repre- 
sentation of the ground and excited state level structure is derived 
from fits to the accurately determined transition frequencies. In 
the first instance, the levels most heavily perturbed by rotational- 
electronic interactions are excluded. Known mass-scaling laws for 
the Dunham coefficients then yield the K coefficients. The non- 
adiabatically perturbed levels were treated in the second instance, 
taking into account wave function mixing of levels in the B'Z+ 
and C'n„ systems. Small adiabatic corrections (at the 1 per cent 
le vel) were included as well. We note here that the values reported 
in lUbachs et af] | |2007|) include additional mass-dependent effects 
on the non-adiabatic interac tion matrix elements that were not in- 
cluded in the initial study of lReinhold et al.l 1 I2OO6I) : differences be- 



Table 1. Catalogue of the most accurate and precise laboratory parame- 
ters for fitting H2 absorptions lines. Represented are all allowed Lyman and 
Werner H2 transitions between the lowest 8 rotational levels in the ground 
and excited states with excited state vibrational quantum numbers up to 20 
and 6 for Lyman and Werner transitions, respectively. The first column pro- 
vides a short-hand notation for the transition: letters denote a Lyman (L) or 
Werner (W) line and the branch, where P, Q and R represent J' - J = -\, 
and 1, respectively, for J and /' the ground state and excited state J- 
levels, respectively; the first integer is the excited state vibrational quantum 
number and the second is J. The second column gives the most precise re- 
ported laboratory wavelengt h and its l-cr uncert ainty and the third c olumn 
provides the reference: 1 = ISaillv et alj )200gh . 2 = lUbachs et al] j2007l) 
(a suffix "a" refers to directly measured wavelengths while "b" refers to 
wavelengths calcu lated from directly m easured lines via combination dif- 
ferences) and 3 = lAbgralletai]h993h for the excited state energy levels 
with ground states derived directly from I Jennings et al] il984) . Note that 
wavelengths with reference 3 are much less precise than those from refer- 
ences 1 and 2. The fourth column gives the oscillator strengths whi c h wer e 
calculated from the Einstein A coefficients given bv lAbgrall eFai] h994l) . 
The fifth column gives the (natural) damping coeffic ients which were cal - 
culated from the total transition probabilities (At) in lAbgrall et al 
The fin al column gives the sensitivity coefficients calculated in lUbachs et alj 
|2007i ) which have estimated uncertainties of typically < 5x 10"'' (see text). 
Only a small excerpt from the full table is presented here; the full table is 
available in the electronic version of this paper 



Trans- A\-jo Ref. f T K 

ition [A] [10"2] [lO's"'] 



LOPl 


1110.062558(3) 


1 


0.05739 


1.87 


-0.0097 


L0P2 


1112.495989(3) 


1 


0.06915 


1.86 


-0.0119 


L0P3 


1115.895530(3) 


1 


0.07381 


1.86 


-0.0149 


L0P4 


1120.248839(3) 


1 


0.07560 


1.85 


-0.0187 


L0P5 


1125.540690(5) 


1 


0.07577 


1.84 


-0.0233 


L0P6 


1131.753504(6) 


1 


0.07481 


1.83 


-0.0286 


L0P7 


1151.6384(14) 


3 


0.07299 


1.82 


-0.0345 



Table 2. Laboratory data for 7 = HD transitions falling in (but not neces- 
sarily detected or fitted in) our spectrum of J2123-0050. The columns have 
the same descriptions as in T able [T]except for the follo wing. T he laboratory 
wavelength references are 4 = lHollenstein et al]j2006l) and 5 = llvanov et alj 
)2008h . The oscillator strengths (/) and damping coefficients (F) were cal- 
culated from the Einstein A coefficients and t otal tr ansition probabilities 
(At), respectively, given bv lAbgrall & Rouefll j200^ . The sensitivity co- 
efficients (K) were calculated bv llvanovetaU 200g) and have estimated 
uncertainties of typically < 1.5 x 10"*. 
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Ref. 
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ition 


[A] 




[10"2] 


[10' s"'] 




LORO 


1105.840555(57) 


4 


0.07436 


1.87 


-0.0065 


LIRO 


1092.001264(58) 


4 


0.29669 


1.76 


-0.0004 


L2R0 


1078.831044(61) 


4 


0.67473 


1.67 


0.0053 


L3R0 


1066.27568(6) 


5 


1.14500 


1.58 


0.0106 


L4R0 


1054.29354(6) 


5 


1.63570 


1.50 


0.0156 


L5R0 


1042.85005(6) 


5 


2.05477 


1.43 


0.0201 


L6R0 


1031.91493(6) 


5 


2.35911 


1.36 


0.0244 


L7R0 


1021.46045(6) 


5 


2.53491 


1.30 


0.0283 


L8R0 


1011.46180(6) 


5 


2.61849 


1.25 


0.0319 


L9R0 


1001.89413(6) 


5 
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1.19 


0.0353 


WORO 
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5 


3.25368 


1.18 


-0.0039 
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Figure 3. Upper panel: The sensitivity coefficients, Kj, for the 7=0-5 HD 
and H2 Lyman and Werner transitions, ;, used in our analysis (dark/coloured 
points) and those not detected or fitted (liglit grey points). The legend ex- 
plains the symbols used in both panels. Lower panel: The distribution of 
transitions with wavelength according to their 7-levels. 



tween the two are nevertheless marginal. We also note that these 
K coefficients are in very good agreement w ith the independent ab 
initio calculations bv lMeshkov et alj ( l2006h : the agreement within 
1 per cent is reassuring for both methods, since it corresponds to the 
estimated uncertainties in both. The K coefficients for the H2 tran- 
sitions are included in Table [T] The K coeffi cients for HD were de- 
rived in ab initio calculations bv llvanov et a l. (2008) and are given 
in Table [2I 

Figure[3](upper panel) shows the K values, which occupy the 
range -0.02 < K < -(-0.03, for all detected molecular transitions. 
Note the general increase in K for bluer Lyman transitions. Adding 
the Werner transitions below A\.jo = 1020 A provides a more com- 
plicated signature of a varying jj. because they shift in the opposite 
direction to the Lyman transitions at similar wavelengths. These as- 
pects of Fig.[3]are particularly important when considering possible 
systematic effects in Section |42l 



3 ANALYSIS AND RESULTS 
3.1 minimization analysis 

If we knew a priori that molecular absorption occurred at a sin- 
gle redshift in a single cloud and no other absorption lines were 
present, measuring A/i/fi would be straightforward: each molecular 
line's redshift could be determined from a fit against a well-defined 
background continuum and its relative deviation from some arbi- 
trary redshift could be plotted against its K coefficient. Equation 
(|TJ implies that A/i/fi would simply be the slope of such a plot. The 
offset would be degenerate with the arbitrary redshift itself (though 
it clearly should be fitted as a free parameter). 

In reality, two main complicating factors are important. 
Firstly, all molecular lines fall in the Lyman-a forest, the series of 



comparatively broad Hi Lyman-o- absorption lines randomly dis- 
tributed in optical depth and redshift bluewards of the quasar's 
Lyman-a emission line. As Figs. [T] and [2] show, many molecular 
lines fall within the absorption profile of one or more forest lines, 
so the latter must also be fitted to provide an effective background. 
Secondly, the molecular lines show 'velocity structure' - several 
absorbing clouds with similar redshifts and different optical depths 
and Doppler widths. The two main spectral features (SFs) high- 
lighted in Fig. |2] are obvious examples, but we demonstrate below 
that each SF comprises more than one absorbing 'velocity compo- 
nent' (VC). 

The 'simultaneous fitting' technique we employ here opti- 
mally and robustly addresses these complications by fitting all de- 
tected Ht and HD absorption profiles, together with the broader 
Lyman-a profiles, simultaneously in a single, comprehensive fit. 
Minimizing x~ allows uncertainties in the forest parameters - the 
effective continuum - to contribute naturally to uncertainties in 
A/j/yU. Other advantages this technique b rings over the earl ier 'line- 
by-line' fitting approach are outlined bv lKing et al.l ( I2OO8I) . 

3.1.1 Free parameters and physical assumptions 

Each absorption line is modelled as a Voigt profile - the convo- 
lution of Gaussian Doppler broadening and the transition's natural 
line-shape. The oscillator strength and natural line- width determine 
the latter; their values are provided in Tables[T]and|2] Three param- 
eters describe each absorption cloud's properties: Doppler width, 
b, column density. A', and redshift, z^bs - This is the only considera- 
tion for the Lyman-a lines but different molecular transitions share 
these parameters in physically important ways. 

For H2 and HD, each /-level has a different ground state popu- 
lation, so all transitions from the same /-level have the same, single 
value of A' in each VC. This is clearly also true for the b and z-jos 
values. However, we make the further assumption that the velocity 
structure is the same in all /-levels, i.e. that a given VC has the same 
Zjibs value in all /-levels. We test the importance of this assumption 
in Section |4~2l and find that it has little impact on our results. We 
also assume that a given VC is characterised by the same, single 
value of b in all /-levels; we test the importance of this assump- 
tion in Section |4n In relating the absorption cloud parameters of 
different transitions in these physically meaningful ways, the num- 
ber o f free parameters is minimized. This is similar to the analy- 
sis in lKing et alj ( |2008|) . though they fitted the molecular oscillator 
strengths as free parameters as well. Finally, because the HD tran- 
sitions are few and relatively weak - recall that only the left-hand 
spectral feature is detected - we further assume that the ratio of an 
HD VC's column density to that of the corresponding VC in the H2 
/=0 transitions is the same for all VCs. As can be seen from Fig.[T] 
there is no strong evidence against this assumption, a conclusion 
corroborated by removing the HD transitions from the analysis al- 
together in Section |4n 

There are several Fen transitions with A\.jo < 1150 A near 
some H2 transitions which, in the z^bs = 2.059 absorber, appear 
quite weak but with velocity structure extending over ~400 km s"' . 
The velocity components in these Feii transitions are marked in 
Fig. [T] This velocity structure was strongly constrained by using 
the Fell /1 1608 A transition, which falls redwards of the Lyman- 
a forest and so has a well-defined continuum. This transition was 
blended with weak, broad Civ /1 1550 absorption associated with 
the quasar itself (Hamann et al., in preparation). The CiV/11550 
absorption was constrained by fitting the corresponding C iv /1 1548 
absorption. Thus, it was necessary to model both metal line species. 
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as shown in Fig.[T] to obtain a reliable model of the Fe n transitions 
appearing near our H2 transitions of interest. The H2 fits are insen- 
sitive to those of the metal lines. 

As noted in Section lTTl a nominal continuum was initially fit- 
ted to the Lyman-Q- forest region before establishing a detailed fit 
to the absorption lines. Simple, low-order polynomial fits connect- 
ing seemingly unabsorbed regions of the spectrum, typically spaced 
apart by >5000 km s"' , provided a qualitatively realistic continuum. 
However, quantitatively different continua were determined by dif- 
ferent authors of this paper; the differences were large enough to 
slightly affect fits to individual absorption lines. And because our 
line fitting approach explicitly connects the properties of molecular 
lines to each other, it is important to recognise that this nominal, 
fixed quasar continuum is unlikely to allow very good fits to the 
many molecular transitions simultaneously. It is therefore impor- 
tant in such cases to fit local continua around transitions where the 
nominal continuum is particularly uncertain. These local continua 
are shown in Figs.[T]and|2l 

The reader will note that, for molecular transitions around 
which broad Lyman-or lines are also fitted, there will be strong 
degeneracies between the Lyman-or line parameters and those of 
the local continuum. For this reason, the local continua have poly- 
nomial degree of, at most, unity (i.e. they are either constants or 
straight lines in wavelength space). Also, the fitting code we use 
indicates cases where the degeneracies are so strong that a local 
continuum fit serves the same purpose as the Lyman-o- line(s). It 
should be emphasised that these ambiguities and uncertainties in 
the local continuum shape all contribute to the uncertainties in the 
molecular line fits in a natural way because, in our^^ minimization 
process, all parameters of the fit are varied simul taneously. Most 
previous works (excluding, e.g., iKing etal.ll2Q08l) efi'ectively de- 
termined fixed, local continua before fitting the molecular lines. 
This means that uncertainties in the continuum fits do not propa- 
gate through to the more interesting parameters of the subsequent 
fit, like A/j/yU. While only a small effect, the uncertainty in A/yZ/j 
would be underestimated with such an approach. 

Weak night sky emission is nominally subtracted during the 
initial flux extraction process, but systematic uncertainties in the 
zero flux level can and do remain. There are many obvious ex- 
amples in Fig. [T] of saturated Lyman-d features which have non- 
zero average flux (either positive or, usually, negative) in their flat- 
bottomed line cores (e.g. features near 3107, 3137, 3171, 3317 and 
3350 A). We therefore included the zero level as a free parameter 
when fitting regions of spectrum which included nearly saturated 
absorption lines (either Lyman-O' and/or molecular). 

We add to our list of free parameters the oscillator strengths 
for the 10 H2 transitions with 7=0-3 falling in the wavelength range 
/lobs = 3345-3415 A. While establishing our fit, these transitions 
had noticeably higher optical depths than predicted by the model. 
For this reason the oscillator strengths for these transitions were left 
as free parameters to be determined in the;^'- minimization process. 
Their fitted-to-calculated oscillator strength ratios, fnt/fcak, deter- 
mined using our fiducial model (see below), range from 1.5 to 2.0. 
Although it is possible that the published oscillator strengths for 
these transitions, which (along with all other molecular lines used 
in our fit) are calculated rather than experimentally measured, are 
incorrect, we consider this very unlikely because these transitions 
are not affected by interference between multiple interacting H2 
states. Fits to H2 transitions, which included some of these partic- 
ular transitions, in other quasar absorbers also suggests their pub- 
lished oscillator strengths are reliable (J. King, private communica- 
tion). An alternative explanation is motivated by the fact that these 



transitions all fall on the combined O vi/Lyman-y8 quasar emission 
line at Aiab ~ 1025-1045 A: perhaps there is some inhomogeneity 
in the H2 column density on the scale of the quasar broad emission- 
line region. Independent evidence for this scenario does not appear 
readily available in our spectrum, so it is difficult to confirm and 
will be explored in future work. Nevertheless, given this possible 
explanation, it is important to test whether these transitions affect 
our final results, and we conduct this consistency check in Section 

ED 

After establishing the best-fitting values of the model param- 
eters described above, we introduce the last free parameter - the 
one of primary interest here - Ayu/yU. Note that only a single value 
of A/j/yU describes all molecular transitions, including all their con- 
stituent VCs. Equation l[T) states that a change in yu manifests itself 
as a pattern of relative shifts between the different molecular transi- 
tions. Note that Fig.|3]shows that the diversity of K values ensures 
little covariance between Ayu/yU and other parameters; the pattern of 
shifts is not degenerate with, for example, the individual VC red- 
shifts. 

3.1.2 Minimizing 

The absorption model is constructed with the free parameters de- 
scribed above and the x" between it and the spectral data is mini- 
mized using the program vpfitQ. This is a non-linear least-squares 
minimization package written specifically for modelling com- 
plex, possibly interrelated absorption lines with a series of Voigt 
profiles. The model fit is convolved with an instrumental resolution 
function for comparison with the real spectrum; we used a Gaus- 
sian function with FWHM = 2.7 km s"' to match the resolution 
inferred from the extracted ThAr exposures. The software 'under- 
stands' that many molecular transitions can arise from the same 
/-level (i.e. that they have the same A') and the ability to link phys- 
ically related parameters, such as b and za^ values, across many 
transitions is inherent in its design. 

t^HlyL was recently added to the vpfit code as a free parameter. 
The operation of a varying yu on the molecular lines is functionally 
the same as the action of a varying a on metal lines. That vpfit re- 
turns the correct values of Aa/a and its uncertainty has been tested 
with a variety of simulations, including models with many strongly 
overlapping VCs which are, in som e cases, blended with transitions 
from unrelated absorption clouds jMurphv et alj|2003l) . However, 
the fit we conduct here is substantially more complicated (though 
operationally similar) to such metal-line fits because the number 
of transitions and links between fitted parameters is much larger. 
Appendix|A]describes a Monte Carlo test of vpfit using our actual 
fit to the molecular, Lyman-o- forest and metal absorption lines in 
J2123-0050 as a synthetic spectrum, to which noise is added and 
into which a value of Ayu/yu = -1-5 x 10"^ is inserted. The^- mini- 
mization of the fit is then re-run on 420 such realizations. The mean 
and distribution of the Ayu//j values recovered agrees well with ex- 
pectations. We also conducted convergence tests using our fit to the 
real data to ensure that vpfit arrived at the same value of Ayu/yU given 
a range of starting values; the default for all the results quoted in 
this paper was zero, but starting values of, e.g., ±10 and ±5 x 10"* 
gave the same results. The self-consistency checks conducted in 
Section|4]also corroborate the reliability of the algorithm. 

The l-cr uncertainties on the best fitting parameters are de- 
rived by VPFIT from the appropriate diagonal terms of the final pa- 

■* http://www.ast.cam.ac.uk/~rfc/vpfit.html. 
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rameter covariance matrix. Given a particular absorption model, 
these errors represent only the formal statistical uncertainties de- 
rived from the flux error arrays of the fitted spectral regions. They 
do not reflect additional sources of error, like wavelength calibra- 
tion uncertainties (see Section |4}. Nor do they represent 'model 
errors' - the possible systematic errors involved with selecting an 
appropriate absorption model. We discuss the selection of the 'best' 
model at length in Section [X2l vpfit assumes that the fluxes (and 
flux errors) in neighbouring spectral pixels are uncorrelated. This 
is certainly not true in practice because the different quasar expo- 
sures were re-binned onto the same final wavelength grid for com- 
bination into a single spectrum. The correlation is reduced by av- 
eraging these exposures. And we are unable to detect the effect by 
comparing the RMS flux variations in unabsorbed regions with the 
final error spectrum; these values are very similar. We demonstrate 
in Section l42l that the binning of the different exposures onto the 
same wavelength grid produces a small systematic effect on A/j/yU. 



3.2 Fiducial absorption model 

Nearly all the H2 and HD transitions detected in our spectrum of 
J2123-0050 were utilized in our fit; some transitions' profiles were 
heavily blended with either nearby Lyman-a lines and/or suspected 
metal lines from one or more absorption systems (e.g. the LOPl Ht 
transition falling at /lobs ~ 3396.2, see Fig.lTJ. No reliable fit could 
be obtained so we chose not to include these molecular lines in our 
sample. 

In total, 86 H2 transitions and 7 HD transitions are included 
in the fiducial fit. Only small regions around the molecular lines 
were fitted: each region contained only enough of the surround- 
ing Lyman-a forest absorption to define the effective continuum 
against which the molecular lines were fitted. This meant that some 
molecular lines were fitted together in the same spectral region, 
joined by fitted Lyman-a forest lines and, in most cases, with a sin- 
gle (straight-line) local continuum and/or zero level. All the fitted 
regions are shown in Fig.[T] Assuming a simple, two VC model 
for the molecular absorption, an initial model of the Lyman-a for- 
est lines, local continua and zero levels was established and refined. 
This same basic structure was used as a starting point in subsequent 
fits to determine the molecular velocity structure. 

Although the absorption profiles in Fig. |2] clearly show two 
spectral features, it is possible, even likely, that the profiles com- 
prise more than two VCs, each of which would be detectable (in 
the absence of the others) given the large number of molecular 
transitions we observe. The velocity structure best representing 
the molecular absorption profile was determined by fitting models 
with increasing numbers of molecular VCs and selecting the one 
with the smallest per degree of freedom, This is one simple, 
standard method for discriminating between different models in x' 
analyses, tho ugh other (ver y similar) 'information criteria' can also 
be used (e.g. lLiddlell2007l) . Note that while (the minimum) it- 
self must always decrease when more free parameters are added to 
a model, xl begins to increase when the additional parameters are 
not statistically justified. Figure |4] (lower panel) shows xl versus 
the number of molecular VCs in the fit. The 4-component model 
has a lower than the 3-component models and so is statistically 
prefened. 5- and 6-component fits were attempted but vpfit found 
the additional components to be statistically unnecessary. That is, 
during the x~ minimization process those components became so 
weak that^- became insensitive to their parameters, so the compo- 
nents were removed. Equivalently, if they had remained in the fit. 
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Figure 4. A,u//i (upper panel) and P^r degree of freedom, xl (lower 
panel), for different velocity structures characterised by the number of fit- 
ted absorption components. The 4-component fit, highliglited with square 
points, has the lowest xl and is therefore the statistically preferred model. 
Two qualitatively different 3-component fits were possible; their results are 
horizontally offset here for clarity. Note the very different values of A/j//i 
they return, exemplifying the inherent systematic eft'ects associated with 
'under-fitting' the absorption profile. The error bars represent l-o" statisti- 
cal uncertainties only; systematic errors are discussed in Section l4!2l 



the final, reduced xl would have been larger than for the fiducial 
4-component fit. 

Figures [T] and |2] show our fiducial 4-component fit. Table |3] 
provides the molecular cloud properties - zjos^ b and the column 
densities in H2 and HD - and their formal 1-cr statistical uncertain- 
ties in all relevant y-levels for each of the 4 VCs. Milutinovic et 
al. (in preparation) and Tumlinson et al. (in preparation) consider 
the absorption cloud properties in detail. 

When visually assessing the fit in Figs. [T] and |2] one clearly 
cannot readily "see" 4 VCs and so one may worry that too many 
components are being fitted. However, this is merely a failing 
of the human eye to recognise coherent structures in almost 100 
transitions simultaneously. For example, it is simple to simulate 
a similar molecular absorption profile as in the real data which 
comprises 4 components, one or more of which might be sub- 
tle enough to go unnoticed when visually assessing several tran- 
sitions but which is statistically significant when ~100 transitions 
are analysed properly. And so we must appeal to the objective, 
well-tested and well-understood model selection technique dis- 
cussed above to recover all the statistically significant structure 
in the line profiles, otherwise we should expect systematic er- 
rors in our results. Previous works have demonstrated that 'under- 
fitting' - fitting too few VCs - causes strong systemati c enws in 
quantities like Aq'/q' l lMurphv. Webb & Flambauml2008l) and Ayu/^ 
jMurphv et al.ll200^ when derived from individual absorption sys- 
tems. 



3.3 Fiducial result 

Given an absorption model for which x^ has been minimized, de- 
termining A/j//j is straightforward: it is simply among the many 



© 2010 RAS, JVDSIRAS 000,[T]{T5] 



Constraint on the variation of jx 9 



Table 3. Molecular absoq^tion line parameters and l-cr statistical uncertainties for the 4-component fiducial fit in Figs.[T]and|2] Note the assumptions discussed 
in Section [3. 1 . 1 I which minimise the number of free parameters. In particular, the ratio of the HD (/ = 0) velocity component column densities to those of the 
corresponding H2 /=0 components is the same for all components. We also provide the total column density for each H2 J level and for HD. 



Com- 




b 






logA'(H2) [cm-2] 






logA'(HD) 


ponent 




[kms-'] 


J = 


/ = 1 


7 = 2 


J = 3 


7 = 4 


7 = 5 


[cm-2] 


1 


2.0593276(5) 


5.14±0.11 


15.06 ± 0.03 


15.38 ±0.06 


15.01 ± 0.03 


15.03 ± 0.02 


13.98 ±0.04 


13.76 ± 0.06 


12.95 ± 0.03 


2 


2.0593290(4) 


1.92 ±0.06 


15.80 ±0.40 


17.52 ± 0.04 


16.23 ±0.16 


15.16 ±0.25 


13.45 ± 0.09 


12.78 ±0.37 


13.69 ± 0.05 


3 


2.0595264(76) 


9.66 ± 0.90 


13.90 ±0.08 


14.30 ± 0.07 


14.10 ±0.05 


14.11 ±0.07 


13.31 ±0.22 


13.47 ±0.15 




4 


2.0595597(8) 


4.02 ±0.16 


13.81 ±0.06 


14.69 ± 0.03 


14.26 ± 0.04 


14.39 ± 0.04 


13.41 ± 0.09 


13.22 ±0.15 




Total 






15.88 ±0.14 


17.53 ± 0.04 


16.26 ±0.11 


15.46 ± 0.04 


14.23 ± 0.03 


14.04 ± 0.04 


13.77 ±0.03 



free parameters determined in the;^'^ minimization process. For the 
4-component fiducial model, we find that 

A^/;/ = (-h5.6±5.5,,„)x 10-^ (2) 

Again, we emphasise that the l-cr uncertainty derives only from 
the spectrum's photon statistics and is calculated from the relevant 
diagonal term of the final parameter covariance matrix. Given the 
range of K coefficients involved (~0.05; see Fig. [3](, equation (TJ 
implies that this uncertainty corresponds to a velocity precision of 
~80 m s"' , or ~0.06 spectral pixels. 

Figure |4] (upper panel) shows the values of A/j//j for different 
models of the molecular velocity structure as characterised by the 
number of fitted absorption components. It graphically illustrates 
the dangers of 'under-fitting'. The two different 3-component fits 
possible in this system give quite different A/i/fi values. This large 
difference should be placed into the context that, since all fits were 
carried out on the same transitions in the same spectrum, the results 
from the different models are by no means independent; the statis- 
tical error bars are not the gauge of "significant" differences here. 
Also, the 2-component model may seem statistically acceptable at 
first since its xl ~ ^ - However, its xl is much larger than for the 4- 
component fit and so it cannot be preferred in any objective sense. 
This is also true of the two different 3-component fits. 

The inadequacies of 2- and 3-component fits are further ex- 
posed in Fig.|5] For the 24 relatively unblended Ht transitions, the 
residuals between the data and the model fit were normalized by the 
flux error arrays, shifted to a common velocity scale and averaged 
together. This forms a 'composite residual spectrum' as a diagnos- 
tic for under-fitting. By combining the residual structure between 
the model and the data for many transitions, we may now better 
appreciate "by eye" how well the model describes the statistically 
significant structure in those absorption profiles. Figure [5] clearly 
shows that the 2-component model fails to reproduce the real ab- 
sorption profile shape, leaving many-pixel excursions outside the 
expected residual range as evidence of additional VCs. This is also 
true for the 3-component models. By comparison, the 4-component 
model leaves no obvious evidence for unmodelled, statistically sig- 
nificant structure. 



4 INTERNAL CONSISTENCY AND SYSTEMATIC 
ERRORS 

4.1 Consistency tests 

The fiducial model in Section \3?2\ includes all detected molecular 
transitions which could be fitted reliably; the data quality and de- 
gree of Lyman-or blending naturally weights each one's contribu- 
tion to Afi/fi. The main advantage of this approach is objectivity: 
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Figure 5. Composite residual spectra (CRS). formed from 24 relatively un- 
blended H2 transitions. Bottom panel: An example H2 transition fitted with 
the 2- and fiducial 4-component models (dashed and solid curves, respec- 
tively). Only a small dift'erence is noticeable by eye. However, the CRS 
panel above shows how significant the additional components are. Upper 
panel: CRS for the 4-component fiducial model (solid line) and the 2- and 
3-component models (dashed and dotted lines, respectively; for clarity only 
one of the two possible 3-component models is shown here). Note the large, 
many-pixel excursions outside the ±1-0" range for the 2- and 3-component 
models, indicating unmodelled velocity structure. No such features exist for 
the 4-component model, as expected. 



we do not select which transitions are 'best' to fit, or those which 
we expect, a priori, to provide the strongest constraints on Aji/p. 
Another approach is to fit only the 'least blended' or the 'well- 
defined' transitions. Deciding which transitions fit these descrip- 
tions is clearly subjective. Nevertheless, as a consistency check, 
we fitted spectral regions containing either no detectable Lyman-a 
blending or very simple and weak blending. We also excluded tran- 
sitions blended with the Fe n lines described in Section l3.1.1l With 
53 H2 and 5 HD remaining transitions, a 4-component model was 
statistically preferred, returning Ayu//j = (+5.0 ± 5.9) x 10"*. 

The 7 = 4 and 5 H2 transitions are weak, typically absorb- 
ing < 30% of the local continuum. It is possible, though we would 
argue unlikely, that unmodelled structure in the nearby Lyman-a 
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forest lines and/or local continua for these transitions will have a 
greater relative effect on the molecular profile fits in these cases. 
A 4-component fit where the 28 7 = 4 and 5 transitions are pre- 
vented from directly influencing A^/yU, effectively by setting their 
K values to zercH, yielded Ayu//i = (+4.8 ± 5.6) x 10"*. The very 
small increase in the statistical uncertainty compared to equation 
^ is symptomatic of how weakly these 7 = 4 and 5 transitions af- 
fect t^iiln, due primarily to their low optical depth. A similar prob- 
lem may exist for the HD transitions, which are also very weak; 
a 4-component model with the HD transitions' dependence on yu 
removed returned A/i/yU = (+4.1 ± 5.8) x 10"*. 

In Section [3.1.1| we described how the predicted optical depths 
for the 7=0-3 H2 transitions falling at A„f,^ = 3345-3415 A 
(i.e. where the Lyman-j8/0 iv quasar emission line falls) were too 
low, possibly because there is some inhomogeneity in the H2 col- 
umn density on the angular scale of the quasar broad line region. 
We therefore allowed these transitions' oscillator strengths to ef- 
fectively vary as free parameters in the fit. This may not be an ac- 
curate parametrisation so we must assess the extent to which it af- 
fects our conclusions. Firstly, if we remove the /j-dependence of 
the 10 7=0-3 Ht transitions concerned (by setting their K values 
to zero), a 4-component fit gives A^/;/ = (+7.0 ± 6.7) x 10"*. An- 
other test is to fit all H2/HD transitions with oscillator strengths 
as free parameters, i.e. the column densities of different transi- 
tions from the same 7-level are allow ed to be different, as in 
the fits conducted bv lKing et al.l l [20o3) . A 4-component fit yields 
A^/// = (+5.4 + 5.5) X 10"*. Notice that the statistical uncertainty 
does not increase discernibly here, even though the number of fit- 
ted parameters has been increased significantly; this is because the 
many column-density parameters for the molecular velocity com- 
ponents of each transition have very little influence on the fitted line 
positions, which is what primarily affects Ayu//j. 

Our fiducial model assumes that a given molecular VC has the 
same fc-parameter in all 7-levels. One can envisage different physi- 
cal effects which may invalidate this model assumption. And while 
we do not find obvious evidence to the contrary in our spectra, it 
may be that our fitting approach - where we make the assumption 
and then build up the fit to the molecular lines by adding VCs until 
all the statistical structure in the molecular line profiles is modelled 
- is not the best method for identifying such an effect. We have 
therefore relaxed the assumption to test its effect on the fitted value 
of Ayu//j: when allowing VCs to have different ^-parameters in dif- 
ferent 7-levels, Ayu/yU = (+7.9 ± 5.6) x 10"*. While there is small 
change in the value of A/j/yU with respect to the fiducial value in 
equation (|2j, it does not represent a strong model dependency. We 
note in passing that allowing VCs to have different /^-parameters in 
different 7-levels does not solve the problem, discussed above, of 
some H2 transitions appearing stronger in the spectrum compared 
to the prediction of the fiducial fit; those transitions come from a 
range of different 7-levels. 

Finally, vpfit allows different values of A/i/yU to be fitted for 
the different H2 7-levels. In this test we maintain the assumption 
that all 7-levels have the same velocity structure, i.e. all 7-levels are 



^ At first it may seem that a simpler, more easily interpretable alternative 
would be to completely remove the 7=4 and 5 (or HD) transitions to test 
their influence on the fiducial result. However, as we describe in Section 
13. 1.1 l and as is apparent in Fig.[T] several molecular transitions are often fit- 
ted together in one contiguous fitting region. This is because different transi- 
tions can be blended together, or be effectively joined by the same Lyman-a 
forest lines. Thus, in principle, removing individual molecular transitions 
from the fit is not possible without removing entire fitting regions. 
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Figure 6. Left panel: t^iilii for each H2 7-level and HD, assuming they all 
have the same 4-component velocity structure. Right panel: t^nlji for two 
groups of transitions, HD plus H2 7=0-1 and H2 J=2-5. The dark/blue 
points derive from a model assuming the same velocity structure in both 
groups, whereas the light/red points (offset to the right for clarity) allow the 
velocity structure for the two groups to differ. Both groups were fitted with 
a 4-component model in both cases. 



fitted with 4 VCs, with each VC having the same redshift in all tran- 
sitions (we relax this assumption in Section |4T2] l. We also derived 
a separate value from the HD lines. The resulting values of A///^ 
are plotted in Fig. [6] (left panel). This plot is useful in clarifying 
the relative contributions the series of transitions from different H2 
7-levels and HD make to the final uncertainty on Ayu/yu. Of course, 
the velocity structure was assumed to be the same in all cases, so 
the different A///yu values are not strictly independent. Nevertheless, 
we note general agreement, especially for the H2 7= 1-3 transitions. 
Since the error bars are relatively large for the other 7-levels and 
HD, we also derived less uncertain values by grouping the HD and 
H2 7=0-1 transitions together. A similar grouping of the H2 7=2-5 
transitions gave a similar value of Ayu/yu. These values are plotted 
on Fig.[6](right panel). 

The above consistency checks indicate that the fiducial value 
of Ayu/yU is not strongly dependent on the assumptions made in our 
4-component fit. Individually relaxing these assumptions produces 
little change in the measured value of Ayu//j. This is also consis- 
tent with there being no important problems with our data, or our 
analysis. 

4.2 Systematic errors 

4.2.1 Specific, known wavelength calibration errors 

Since a varying ^ should cause velocity shifts between molecular 
transitions, wavelength calibration errors which distort the wave- 
length scale are clearly important potential sources of systematic 
errors in Afi/ii. Note that, as implied by equation l|T), effects which 
shift the wavelength scale by a constant velocity at all wavelengths 
are generally unimportant. The first, most obvious possibility is 



This is not strictly true when many quasar exposures are combined, as is 
the case here. If different velocity shifts are applied to the different quasar 
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that the wavelength scale of the ThAr exposures was incorrectly de- 
termined. This is easily checked. The wavelength calibration resid- 
uals, with RMS ~ 80 ms"', are consistent with being symmetri- 
cally distributed around the final wavelength solution at all wave- 
lengths. Because many H2/HD transitions are used, with a range of 
wavelengths, this reduces the overall calibration error substantially. 
For example, over the 12 HIRES echelle orders covering the wave- 
length range where the molecular lines fall (3071-3421 A), 150- 
200 ThAr lines are used to calibrate our ThAr exposures. Clearly 
though, systematic trends in these r esiduals are sti l l poss ible. Us- 
ing the same technique employed in iMurphv et all ( |2007|) to track 
systematic patterns in the ThAr calibration residuals, the possi- 
ble distortion in the wavelength scale between 3070 and 3430 A 
is <30ms"'. This corresponds to a systematic error in of 
±2.0 X 10"* almost. 

Once determined from the ThAr exposures, the wavelength 
solutions are simply applied to the corresponding quasar exposures. 
Drifts in the refractive index of air inside HIRES between the ThAr 
and quasar exposures will therefore cause miscalibrations. The 
temperature and atmospheric pressure drifts during our o bserva- 
tions w ere < 1 K and < 1 mbar respectively. According to the lEdlerJ 
( Il966h formula for the refractive index of air, this would cause dif- 
ferential velocity shifts between 3070 and 3430 A of <10ms"', 
which are negligible. 

The quasar and ThAr light traverse similar but not identical 
paths through HIRES. Only the quasar light passes from the tele- 
scope into HIRES and the ThAr light (nearly) uniformly illumi- 
nates the slit while the quasar light is centrally concentrated. Differ- 
ent distortions of the wavelength scale - within individual echelle 
orders ('intra-order') and over longer wavelength ranges ('long- 
range') - may therefore occur in the quasar and corresponding 
ThA r exposures. 

iGriest et"al] l l201Clh recently identified intra-order distortions 
in HIRES. They compared the wavelength scales established us- 
ing ThAr exposures with those imprinted on quasar exposures by 
an iodine absorption cell. The distortions are such that, for differ- 
ent transitions at the same redshift in a quasar spectrum, those at 
the order edges appear at positive velocities with respect to tran- 
sitions at the order centres when calibrated with a ThAr exposure. 
The peak-to-peak velocity distortion is ~500ms"' at ~5600 A and 
may worsen ( improve) with inc reasing (decreasing) wavelength; 
see figure 4 of lGriest et al] j201(]|) . The effect is similar (though not 
identical) for all echelle orders in the wavelength range ~5000- 
6200 A covered by the iodine cell absorption. If we assume that 
similar intra-order distortions apply to our spectra of J2 123-0050 
at much bluer wavelengths then, because the molecular transitions 
of interest lie at different positions along different echelle orders, 
we should expect the effect on A/j/yU to be suppressed. To illus- 
trate this, a very crude estimate of the 'residual' velocity distortion 
is just the observed ~500ms"' peak-to-peak intra-order value re- 
duced according to the number of molecular transitions observed, 
i.e. ~ 500ms"V V93 » 52 ms"'. This corresponds to a systematic 
error in A/j/yU of approximately ±3.5 x 10"*. 

This crude calculation fails to take into account the strong 
variation in the spectral SNR across the wavelength range contain- 
ing the Hi/HD transitions, their different sensitivities to yu variation, 
where they fall with respect to echelle order edges in our particu- 



exposures, and if the relative weights of the exposures vary with wavelength 
when forming the final, combined spectrum, then small relative velocity 
shifts will be measured between transitions at different wavelengths. 



lar spectrum and other factors such as the degree of blending with 
Lyman-a forest lines of each molecular transition. A simple Monte 
Carlo simulation which takes into account the first two of these 
effects - SNR and K for each transition - was undertaken as fol- 
lows. Each realisation comprised a velocity shift, Av,, chosen ran- 
domly for each transition, i, from the 500 m s"' peak-t o-peak inter- 
val sp anned by the intra-order distortions identified bv lOriest et al.l 
l l20ia) . The square of the SNR of the spectral data surrounding the 
molecular transitions then weighted a linear least squares fit of the 
Av, /c versus K values, the slope of which provided a value of A;///j 
for that realisation - see equation (TJ. The distribution of values for 
Ayu//j over hundreds of thousands of Monte Carlo realizations was 
close to Gaussian with an RMS of 4.8 x 10"*. 

However, the most direct calculation of the possible system- 
atic effect on Ayu//j, which takes into account all remaining effects 
(e.g. positions of transitions with respect to order edges, line blend- 
ing etc.), can be obtained by attempting to remove the intra-order 
distortions from our individual quasar exposures. To this end we 
applied a -250 m s"' velocity shift to the wavelengths of all echelle 
order centres for all our quasar exposures. The shift was reduced 
linearly with distance from the order centres to reach -1-250 m s"' at 
the order edges. The spectra were recombined again to form a final 
1 -dimensional spectrum to which we fit the 4-component model as 
before. The value of A;u//j derived from this 'corrected' spectrum 
was (+3.7 ± 5.5) x 10"*. Compared with the fiducial 4-component 
result of (+5.6 ± 5.3) x 10"*, this represents a smaller effect, ap- 
proximately ±1.9 X 10"*, than expec ted from the cruder estimates 
above. Given that lGriest et alj ( l201(j) observed smaller distortions 
at shorter wavelengths, it may also overestimate the real effect. 

To summarise the above considerations, i f the intra-order dis- 
tortions of the wavelength scale identified by iGriest et"aL 
at 5000-6500 A also affect the much bluer wavelengths of inter- 
est here, the implied systematic error on Ayu/^ is approximately 
±1.9 X 10"* in our particular spectrum. However, the main ass ump- 
tion in this estimate is that we have modelled the lGriest et akl intra- 
order distortions reliably. The physical explanation for those distor- 
tions is not yet known, so we do not know whether, for example, the 
phase of the saw-tooth distortion pattern with respect to the order 
edges is the same for all quasar exposures taken at different times 
and with different telescope and observational conditions (e.g. tele- 
scope pointing direction, seeing etc.), nor do we know whether the 
intra-order distortions follow a general saw-tooth pattern at all, or 
if the distortions in the UV can be predicted based on those in the 
iodine cell's calibration range. These important questions should 
be resolved with future observations and careful re-calibrations of 
HIRES. Thus, our estimate of ±1.9 x 10"* for the systematic error 
due to intra-order distortions may also be subject to model errors. 

4.2.2 General, unknown wavelength calibration errors 

In Fig.[3]one observes a general decrease in the K coefficients for 
the fitted Lyman H2 transitions with increasing wavelength. A long 
range, monotonic distortion of the wavelength scale is therefore 
nearly degenerate (though not completely) with the effect of a shift 
in yu on the Lyman lines. Put another way, a systematic effect which 
produced a long-range wavelength distortion would cause a strong 
systematic effect on Ayu/yu if only the Lyman H2 transitions were 
fitted. It is therefore important to fit the Werner transitions as well: 
their K values are very different to those of the Lyman transitions at 
similar wavelengths, potentially breaking the degeneracy between 
Ayu/yU and long-range wavelength distortions and improving resis- 
tance to systematic errors. If we remove the yu-dependence of the 
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12 H2 and 1 HD Werner transitions in our fit, a 4-component model 
is statistically preferred and returns A////J = (+5.2 ± 5.6) x 10"*, 
very similar to the fiducial result in equation l|2j. 

However, the SNR at bluer wavelengths where the Werner 
transitions fall is much lower than for the reddest Lyman transi- 
tions, so the effect of the Werner transitions in breaking this de- 
generacy between A/j//j and long-range wavelength distortions is 
reduced. To illustrate this we inserted such a distortion into a sim- 
ulated version of our spectrum. The wavelength scale for one of 
the Monte Carlo realizations from Appendix A was compressed 
according to 



Avj 



(3) 



where Aj is the initial (observed) wavelength of a given pixel, Avy 
is the velocity shift applied to it, (A/j/yu)sys is the systematic error 
in Ayu/// we seek to mimic, and the constants a and b are set to 
mimic the general decrease in K with increasing wavelength for 
the Lyman transitions (see Fig. [3]l. To obtain a spurious shift in 
A/J//J of approximately +15 x 10"'', we set (A/j//j)sys = 15 x 10"*, 
a = -1.3 X 10"* A"' and b = 0.4345. Indeed, when this spectrum is 
fitted with the fiducial model used to generate it, Afi/ji shifts from 
4.0 X 10"* before compression to 16.2 x 10"* after compression, 
as expected. And when the //-dependence of the 12 H2 and 1 HD 
Werner transitions is removed when fitting the compressed spec- 
trum, Afi/jj only changes by small amount to 17.3 x 10"*. That is, 
simply removing the Werner transitions' dependence on /j is not 
a very effective test for long-range distortions of the wavelength 
scale. 

A more meaningful test is to fit only the Lyman and Werner 
transitions in the bluest part of the spectrum and not the redder 
Lyman transitions. If we fit the 51 Werner and Lyman lines blue- 
wards of /lobs = 3230 A in our simulated spectrum, Afi/fi only shifts 
from 6.1 to 7.4 x 10"* when the compression is introduced, indi- 
cating the expected resistance to this systematic error. With the /i- 
dependence of the 12 H2 and 1 HD Werner transitions removed in 
this fit, A/i/fi is 7.6 and 1 1.3x10"* before and after compression, re- 
spectively. That is, if there is a long-range, monotonic distortion of 
the wavelength scale, we expect to find a substantial shift in Ayu//j 
when the ^-dependence of the Werner transitions is removed. For 
the real spectrum, the value of A/i/yU derived from the Lyman and 
Werner lines bluewards of A^^,, = 3230 A is (+12.1 ± 11.7) x 10"* 
and, when the /i-dependence of the Werner transitions is removed, 
this changes to (+12.8 ± 14.3) x 10"*. Thus, despite this being a 
reasonable test for long-range distortions of the wavelength scale, 
it fails to reveal evidence for such an efi'ect in our spectrum of 
J2123-0050. 



4.2.3 Errors from spectral re-dispersion 



As discussed in Section [2n the spectrum to which we fit our ab- 
sorption profile models is the weighted mean of several individ- 
ual quasar exposures. In combining these spectra, they were all re- 
dispersed onto the same wavelength grid. Note that the data re- 
duction procedures we employed meant that each exposure was 
only re-binned in this way just once. Nevertheless, this rebinning 
introduces correlations between the flux (and flux uncertainty) of 
neighbouring pixels. In principle, this should have a small effect on 
the fitted centroid of any spectral feature and so, even though we 
fit many molecular lines simultaneously, a residual effect on Ap/p 
may exist. To test the possible size of this effect we re-combined 
the quasar exposures onto a variety of slightly different wavelength 



grids with different dispersions (i.e. kms"' per pixel) and wave- 
length zero points. The values of Ap/p derived from these different 
versions of the final combined spectrum varied by ±0.8 x 10"* at 
most. 



4.2.4 Velocity structure 

The fiducial absorption model assumes that the molecular velocity 
structure is the same in all i-leve ls (see Section|3.1.U . In a G alactic 
line of sight with H2 absorption, Ijenkins & Peimbera jl997h found 
that (part of) the absorption profile becomes broader and systemat- 
ically shifts in velocity with increasing J. The absorption profiles 
comprise many velocity components spanning velocities greater 
than the broadening and shifting observed. This is consistent with 
the velocity structure being the same in all y-levels where some 
constituent VCs have low optical depths in the low- J transitions but 
are relatively strong in the higher-7 transitions. Thus, if one knew 
the velocity structure - or one used as many VCs as required to 
model all the statistically significant structure in the observed line 
profiles - one may not necessarily find that individual VCs broaden 
or shift with increasing J. Instead, one would simply find that the 
relative column densities of neighbouring VCs changed with in- 
creasing J. However, if any significant unmodelled structure re- 
mains, and its effective optical depth varies as a function of J, it 
may cause a systematic effect in Ap/p. An exam ple of this effect in 
an H2-bearing quasar absorber was identified bv lNoterdaeme et al] 
( l2007i) . 



Ivanchik et al] ( l2005h considered this possibility, pointing out 



that the potential effect on Ap/fi is diminished because transitions 
of different J are interspersed in wavelength space. This is also ev- 
ident in Fig. [3] (lower panel). Although we emphasise that Figs.|4] 
and |5] demonstrate that all the statistically significant structure in 
the molecular transitions has been modelled, we can make a direct 
test for this effect by allowing Ap/p and the velocity structure to be 
different for different H2 7-levels and HD. In Section|4T|we made 
only the former assumption and found that a stringent test, with 
relatively small statistical uncertainties, was only possible when 
grouping together different H2 y-levels and HD. By introducing 
more free parameters in this new test - the different VC redshifts 
- for different H2 /-levels and HD (or groups thereof), we expect 
even larger uncertainties. So we again form groups of H2 /-levels 
and HD, deriving a single value of Ap/p with a 4-component ve- 
locity structure in the HD and H2 /=0-l levels and another value 
with a different 4-component model for the H2 7=2-5 levels. The 
two values are shown in Fig.[6](right panel): we see no strong evi- 
dence for the effects of unmodelled velocity structure as a function 
of species (H2 /-level or HD). 



4.2.5 Summary of systematic errors 

To summarise the above discussion, by far the two most important 
potential systematic effects are small distortions in the ThAr wave- 
length scale and infra-order distortions of the quasar wavelength 
scale relative to the ThAr exposures. These contribute systematic 
errors in Ap/p of approximately ±2.0 and ±1.9x 10"*, respectively. 
Effects from re-dispersion of the spectra may contribute as much as 
±0.8 X 10"*. 
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5 SUMMARY AND DISCUSSION 

Our final result for the z^^s = 2.059 absorber observed in the Keck 
spectrum of J2 123-0050 is 

A^/M = (+5.6 ± 5.5„„ ± 2. V) X 10-« . (4) 

This includes the l-o" statistical error [see equation ^] and the 
quadrature addition of the three main potential systematic errors 
discussed in Section l42l long- and short-range wavelength cali- 
bration errors and effects from re-dispersion of the spectra. It is 
difficult to estimate the confidence level represented by the quoted 
systematic error component, though it is likely to represent a higher 
confidence (i.e. >68 per cent) than the 1-cr level represented by the 
quoted statistical uncertainty. However, we emphasise that the sys- 
tematic error estimate for short-range calibration errors is based 
only on the limited informatio n about intra-order wavelength dis- 
tortions currently accessible jGriest et all l2010h : it is therefore 
model-dependent and may under- or over-estimate the true effect, 
if present. 

Our result in equation ^ is based on fitting an absorption pro- 
file model to particular regions of the spectrum which contain 86 H2 
and 7 HD transitions, most of which are blended to varying extents 
with broader Lyman-a forest lines, and a few of which blend with 
well-constrained, weak metal absorption. All these lines, together 
with local continua and zero flux level adjustments for most lines, 
are fitted simultaneously to determine a single, best-fit value of 
Ayu//j. To reduce the number of free parameters in the fit, two main 
physically motivated assumptions were made about the molecular 
transitions: that the velocity structure is the same for all y-levels 
and that individual velocity components have the same Doppler 
broadening parameter in all y-levels. Also, since the HD lines are 
so weak, it was convenient to assume that the HD (J = 0) velocity 
components follow the same relative optical depth pattern as the 
H2 y = transitions. Finally, the oscillator strengths of some H2 
transitions were effectively treated as free parameters because their 
optical depths are higher than expected based on our fiducial model; 
see Section [3. 1.11 However, we have individually relaxed each of 
these assumptions and/or removed the transitions they pertain to, 
and find the fiducial result in equation (|4j to be quite insensitive to 
them. 

Figure[7]compares the contemporary astronomical A////J con- 
straints outside our Galaxy. Our new constraint is the first from 
Keck with precision comparable to those from the VLT. Agreement 
is seen between ou r Keck result and the three VLT constraints from 
iKing et al] j2008l) who used similar calibration and fitting tech- 
niques to those described here. Such consistency is encouraging, 
allowing us to (at least formally) derive a weighted mear0 value 
of A^/fi = (+3.5 ± 2.8) X IQ-^ for redshifts z = 2.0-3.1. Statisti- 
cal errors are shown as thinner error bars with shorter terminators 
in Fig.|7] Although the statistical err or in our new mea surement is 
slightly smaller than for two of the iKing et all ( l2008h absorbers, 
it is slightly worse than for the z^bs = 2.811 absorber towards 
Q 0528-250, i.e. 3.9 x 10"^. The main reason for this is the lower 
SNR of our current Keck spectrum of J2123-0050, despite the fact 
that its spectral resolution is higher and more molecular transitions 
were utilized in our analysis. Since J2 123-0050 is relatively bright, 
improvements in both the statistical and systematic uncertainties 
can be achieved with new observations of moderate duration. And 



' We added the statistical and systematic error components on our new 
result in quadra ture for this calculation. No systematic error estimates are 
available for the lKing et alj results. 
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Figure 7. Current extragalactic Afi/fi constraints. The legend gives the ref- 
erence for each point and the method used, i.e. H2 or NH3 quasar absorp- 
tion lines, l-o" statistical uncertainties are indicated with thinner error bai's 
with shorter terminators while the thicker eiTor bars with longer termina- 
tors show the systematic error added in quadrature with the statistical ones. 
Diamonds represent the VLT constraints; our new measurement is the only 
one from Keck. Note that two of the : King et al.| 1 1200^ absorbers use the 
sa me raw VLT data as tho se analysed previously bv llvanchik et al. ( 200^ 
and lReinhold et al.l i2006l) . We omit the results of those two studies assum- 
ing that the improved analysis in iKing et al] has yielded more reliable re- 
sults. For cla iity, we also omit the othe r rec ent re-analyses of the sa me two 
absorbers by I Wendt & Reimerj fcOOSh and [Thompson et al J J2009li : these 
w orks also yi elded null results with somewhat larger errors than reported 
bv lKing et al.L 



since it is located near the celestial equator, J2123-0050 can easily 
be observed from both Keck and VLT, enabling a further opportu- 
nity to check for instrumental systematic errors. We have recently 
undertaken such observations and a future paper will report the re- 
sults. 

At z < 1, comparison of the radio inversion transi- 
tions of NH3 - which have enhanced sensitivity to ^-variation 
jvan Veldhoven et al.l20"04l ; lFlambaum & Kozl(wll2007h - with less 
sensitive molecular rotational lines (e.g. HCO^, HCN), has yielded 
two very strong constraints, A/j/yU = (+0.74±0.47s,;„±0.76sys)xl0"'' 
at z = 0.685 (Murphy et al. 20oi) and (+0.08 ± 0.47,^^) x 10"'* at 
z = 0.889 l iHenkel et al.ll2009l) . These results are plotted in Fig.|7] 
for comparison with the higher redshift H2/HD constraints. It is 
clear that the radio constraints have superior precision and, by cur- 
rent estimates, smaller potential systematic errors. However, direct 
comparison of the radio and optical constraints is difficult because 
of the possibility, in principle, for spatial variations in /i: the dif- 
ferent molecular species (NH3 and H2/HD) trace regions of dif- 
ferent densities and, therefore, different spatial scales and environ- 
ment. If fi does vary, we do not know what that variation depends 
on, so it is presumptuous to prefer one ty pe of measurement over 
the other. Indeed, this is highlighted by iLevshakov et al.l 1 I2OO8I) 
who studied NH3 inversion emission lines from numerous Galac- 
tic molecular clouds. With statistical errors from previous literature 
of ~0. Ikms"' they find velocity offsets between the NH3 inver- 
sion and rotational molecular emission of up to |Av| ~ 0.5 km s"' 
in individual systems. This might indicate spatial variations in fi 
throughout our Galaxy, although intrinsic shifts between emission 
lines of different molecules are to be expected. The possibility for 
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both space- and time-variations in yu is even more important given 
the different redshift ranges currently probed by the radio and opti- 
cal constraints. 

Several other quasar absorption- and emission-line techniques 
for constraining valuations in combinations of fundamental con- 
stants involving are also of note. For example, comparison of 
Hi 21 -cm absorption with metal absorption lines constrains the 
quantity X = gpO-^/yj where gp is the proton g-factor. A variety 
of analyses using different metal-ions hav e been performed (e.g . 
IWolfe et al.l 1 19761 : iTzanavaris et"ai] l2005l : iKanekar et all 12000) . 
with m easurements from 9 abs orption systems at 0.23 < Zabs < 
2.35 by ITzanavaris et al] ( l2007h providing a constraint of AX/X = 
(+6.3 ± 9.9) X 10"^ However, the fact that the H i 21 -cm and metal- 
line velocity structures are observationally dissimilar - probably 
because the radio morphology of most background quasars is not 
point-like - means that improvements must come by averaging 
over many sight-lines and/or carefully selecting quasars with point- 
like radio morphologies. Comparing Hi 21 -cm with the 'main' 
OH 18-cm absorption lines, which constrains F = gpQia^Y-^^ 
jChengalur & Kanekadl2003l) . suffers less from thi s problem be- 
cause of the similarity of the wavelengths concerned. iKanekar et al.l 
( l2005h analysed two absorbers at Zabs = 0.685 and 0.765 to obtain 
the constraint AF/F = (-1-4.4 + 3.6,„a, ± W,y,) x 10"^ Finally, a 
promising technique is to use the 'satellite' OH 18-cm lines which, 
when the ground state level populations are inverted, results in the 
two transitions being 'conjugate' - one occurs in emission, the 
other in absorption, but they have the same optical depth profile 
and are guaranteed to arise in the same gas clouds. Comparing 
the sum and difference of the optical depth profiles of the 1612 
and 1720 MHz conjugate sat ellite lines constrains G = gp(//Q'^)'"'*^ 
dChengalur & Kanekarllioosh . With only two such systems known 
outside our Galaxy, and with their published spectra hav ing very 
low SNR, no definitive constraints on AG /G yet exist jParlinel 
|2004 IKanekar et al]|2005l) . iKanekail ( |2008|) has presented a pre- 
liminary constraint of \AG/G\ < 1 1 x 10"' for a system at z = 0.247 
and also < 12 x 10"' for a Galactic system, Centaurus A. 

While these constraints on variations in X, F and G are im- 
portant in their own right, it is difficult to compare them with con- 
straints on Afi/fi. As we mentioned in Section [Tj predictions for 
the relationships b etween variations in, e.g., yu and a are strongly 
model-dependent jDent et af I[200i). Currently, there is also no 
overlap with the redshift range occupied by the Hi/HD constraints 
on Aji/fi and variations in X, F and/or G. 

A striking feature of Fig. |7| is the paucity of direct yu mea- 
surements outside our Galaxy. Compare this with, for example, the 
143 constraints on a-variation available from Ke ck/HIRES over 
the redshift range z = 0.2-A.2 jMurphv et"aLll2004h . Clearly, much 
larger samples of both NH3 and H2/HD constraints, in overlapping 
redshift ranges, would allow additional tests for systematic errors 
and stronger conclusions to be drawn. In this respect it is impor- 
tant to reduce the statistical errors in individual H2/HD measure- 
ments. This can only be done by substantially increasing the SNR 
of the optical spectra. We have demonstrated here that, at least for 
the Keck/HIRES spectrum of J2 123-0050, the statistical error still 
dominates over systematic errors [equation (|4]l]. Nevertheless, if 
such systematic errors are not also reduced by improving the spec- 
tral SNR, it becomes more important to greatly increase the number 
of H2/HD absorbers for measuring Afi/ji in order to approach the 
precision demonstrated in the NH3 measurements. Surveys target- 
ing particularly gas- and/or metal-rich absorber s have been success- 
ful in discovering more H2/HD systems (e.g . iLedoux et al.ll2003l ; 
IPetitiean et alj2006l ; lNoterdaeme et al.l2008n . The large number of 



damped Lyman-o- systems identifi ed in the Sloan Digital Sky Sur- 
vey (e.g. lProchaska & Wolfel2009l) offers a new avenue for discov- 
ery, J2123-0050 being one such example. Of course, finding more 
NH3 absorbers, particularly at z > 1 would be highly desirable as 
well. 
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Figure Al. 420 Monte Carlo simulations of the fiducial 4 component ab- 
sorption model used in our analysis. The value of A/i//j returned from the 
simulated spectra (vertical dotted line) corresponds well to the input value, 
+5 X 10"^. The mean l-o" error and the standard deviation are shown to be 
consistent. 



duced from the fiducial fit we established to the real data. That is, 
each simulated spectrum is produced simply by taking the fit with 4 
molecular velocity components, plus all the fitted Lyman-o- forest 
lines, interloping metal lines, local continuum and zero flux level 
adjustments, and adding Gaussian noise with a cr of 0.8 times the 
error array of the real spectrum. The artificially high SNR in the 
simulations was used to ensure that Lyman-a blends that were only 
marginally statistically required in the fit to the real spectrum were 
not removed statistically when fitting the simulated versions. 

Each simulated spectrum was then fitted with vpfit with ini- 
tial guess parameters set to be those used to produce the simu- 
lated spectrum. An input value of Ayu/yU = -1-5 x 10"'' was used 
to test whether vpfit recovered that value when started from an ini- 
tial guess of zero. Figure [ATI shows the results of fitting 420 simu- 
lated spectra. The mean l-o" uncertainty derived on individual Ayu//j 
measurements is 4.4 x 10"* which is the expected value given the 
value of 5.5x 10"* from the real spectrum and a scaling of 0.8 times 
its flux errors. This corresponds well to the RMS of 4.1 x 10"^ for 
the 420 Ayu/yu measurements; it may even be that the individual un- 
certainties returned by vpfit are slightly conservative. The mean 
Ayu/yU value retrieved is 4.8 x lO"*. 

This paper has been typeset from a Tj3C/LiTjjX file prepared by the author 



APPENDIX A: MONTE CARLO TEST OF FITTING 
ALGORITHM TO DETERMINE A///^ 

The fiducial fit performed on the spectrum of J2123-0050 to de- 
termine A/j/yU is a fairly large and complex one. It contains more 
fitted spectral pixels and more fitted parameters than most, if not 
all, previous fits of its kind in the literature. Furthermore, the fitted 
parameters are not all free and independent, but are tied together 
in a variety of physically meaningful ways. Therefore, to ensure 
that the fitting code vpfit determines A/j/yU correctly in this con- 
text, we performed a Monte Carlo test using simulated spectra pro- 



APPENDIX B: SUPPORTING INFORMATION 

Figure IBll is the complete version of Fig. [T] It will appear in the 
online version of this paper, not in the printed version. Table Bl, 
the complete version of Table 1, is provided in the online version 
only as a machine-readable ASCII file. 
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Figure Bl. All regions of the J2123-0050 Keck spectrum fitted simultaneously in our analysis. The spectrum (black histogram) is normalized by a nominal 
continuum (upper dotted line) fitted over large spectral scales. Local Unear continua (upper dashed lines) and zero levels (lower dashed lines) are fitted 
simultaneously with the H2/HD and broader Lyman-a lines. The fits are shown with solid grey/green hnes. Hi/HD transitions are labelled and their constituent 
velocity components are indicated by grey/green tick-marks immediately above the spectrum. Higher above the spectrum are tick-marks indicating the positions 
of Lyman-Q" lines (blue) and Fe n lines (red). Note that the metal-line velocity structure is constrained with the Fe n /n608 A transition shown in the final panel 
of the figure. The residual spectrum (i.e. [data] - [fit]), normalized to the 1-cr errors (faint, horizontal solid lines), is shown above the tick-marks. 
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